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Abstract 

We suggest that exotic sphere fibrations can be mapped to band topologies in condensed matter 
systems. These fibrations can correspond to geometric phases of two double bands or state vector 
bases with second Chern numbers m + n and — n respectively. They can be related to topological 
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c/3 I insulators, magneto-electric effects, and photonic crystals with special edge states. We also consider 



time-reversal symmetry breaking perturbations of topological insulator, and heterostructures of 
topological insulators with normal insulators and with superconductors. We consider periodic 
TI/NI/TI/NI' heterostuctures, and periodic TI/SC/TI/SC heterostuctures. They also give rise to 
models of Weyl semimetals which have thermal and electrical transports. 



I. INTRODUCTION 

In this article, we propose to realize the exotic spheres as the geometric phases in con- 
densed matter systems. We suggest that the exotic sphere fibration can be realized as the 
geometric phases in condensed matter systems, cold atomic or molecular systems. 

The topological description of the quantum states of matter gives a new method in 
describing the condensed matter. Condensed matter systems with band structures that 
have nontrivial topological properties give a new type of materials with properties that are 
robust under certain perturbations. There are many interesting topological properties for 
example the appearance of edge states and the existence of gapless surface states. 

The topological insulator has a bulk gap, while it has topologically protected edge states 
on the boundary of the topological insulator. Thus, the topological insulator is an insulator 
in the bulk while has gapless edge states on the boundary of the topological insulator. 
The 3D topological insulators have been shown in materials for example,-"- Bi2Te3, Bi2Se3, 
Sb2Te3, Bii„a.Sbj:. The 2D topological insulators have been observed in HgTe quantum 
wells. -1^ 

For the 3D topological insulators, the topologically protected surface state realizes itself 
by the non-trivial spin texture on the surface band of the topological insulator. The topo- 
logical surface state is protected by time-reversal symmetry. The existence of odd number 
of or single surface Dirac point is robust in the presence of nonmagnetic impurities, and 
other time-reversal symmetry preserving perturbations. The difference between topological 
insulator and normal insulator can be distinguished by a Z2 invariant.-"- The 3D topological 
insulators can be described^ by topological field theory with a 9 variable. Due to the time 
reversal symmetry, 9 takes values of or tt, modulo 27r. The 9 term gives magnetoelectric ef- 
fects, with magnetoelectric polarization in the materials.— lii The spin polarization of surface 
band and magnetoelectric polarization of the materials, can be experimentally measured. 

In early days of differential topology, John Milnor constructed a seven dimensional com- 
pact space called exotic sphere. The space he constructed has the property that it has a 
continuous one to one map to the round sphere, and yet it cannot be mapped to the round 
sphere smoothly. There are 28 Milnor exotic spheres that are mutually distinct from each 
other. The construction was based on distinct bundles over the four dimensional manifold. 
Bundles over four dimensional manifold may be used to describe geometric phases in con- 



densed matter and atomic systems. It is interesting to see whether one can build the Milnor 
exotic spheres into the theory of condensed matter and atomic systems. It is also nice to 
connect subjects in mathematics to subjects in theoretical physics. 

In this paper, we suggest that exotic sphere fibrations can be mapped to band topologies 
in condensed matter systems. These fibrations can be mapped to geometric phases of two 
double bands or state vector bases with their associated second Chern numbers. They can be 
related to topological insulators, magneto-electric effects, and photonic crystals with special 
edge states, among other aspects. It is nice to understand the physics of topological insulator 
in the situations when it is placed adjacent to other kinds of materials. We also consider 
time-reversal symmetry breaking perturbations of topological insulator, and heterostructures 
of topological insulators with normal insulators and with superconductors. 

The organization of this article is as follows. In section |Tll we discuss that exotic sphere 
fibrations can be related to band topologies in condensed matter systems. These fibrations 
can be mapped to geometric phases of two double bands or state vector bases with the 
second Chern numbers m + n and —n respectively. In section IIIIl we discuss their rela- 
tion to topological insulators, and magneto-electric effects. In section HVl we also discuss 
their relation to photonic crystals with special edge states. In section |Vl we consider time- 
reversal symmetry breaking perturbations of topological insulator, and heterostructures of 
topological insulators with normal insulators, and periodic TI/NI/TI/NI' heterostuctures. 
In section IVll we consider heterostructures of topological insulators with superconductors, 
and periodic TI/SC/TI/SC heterostuctures. These structures also give rise to models of 
Weyl semimetals which have thermal and electrical transports. In section IVIIl we also sug- 
gest relevance to other possible materials such as cold atom systems and semiconductor 
systems. 

II. GEOMETRIC PHASES OF TWO DOUBLE BANDS OR STATE VECTOR 
BASES AND RELATED TOPOLOGIES 

One of the interesting types of fiber bundles are 3-sphere bundles over 4-sphere. Such 
fibrations can be constructed by patching two -R^ x S^ and identify their overlapping region 
by a diffeomorphism. One can divide 5*^ into three regions: a north patch Rfi)', a middle 
patch [— e, e] x S^; and a south patch R^)- The Rf^-. x S^ can be parametrized by a quaternion 



u and a unit norm quaternion v, in which the u belongs to the i^ts and the v belongs to 
the 5*^. Similarly, the i?!) x S^ can be parametrized by a quaternion u' and a unit norm 
quaternion v'. The transition function is defined on the middle patch, and it is 

||m|| denotes the norm of the quaternion u, while \\v'\\ =\\v\\ = 1. Such fibrations can be 
classified by two integers (n + m, —n). 

The fibration of S^ over S'^ can be characterized by the map from the middle patch S^ 
to the structure group 5*0 (4) which corresponds to the rotational symmetry of the S^ fiber. 
This map is characterized by homotopy group 713(5*0(4)) = Z © Z, where 7!'3{SU{2)) = Z, 
and so(4) = su(2)(i) x sn(2)(2). This fibration is characterized by two integers, which 
correspond to {n + m, —n). The two integers n + m and —n correspond to the second Chern 
numbers of the s-u(2)(i) and su{2)(2). This may be viewed as 

4'^ = At /" (trF« A F«) = n + m, (2) 

(^) = J_ f (trF(2)AF(2)) = -n, (3) 



and the integration is on the four dimensional base manifold. 

This can be interpreted a.s n + m instantons of the sm(2)(i) gauge fields, and n anti- 
instantons of the 5^(2) (2) gauge fields. This can also be viewed as m instantons of s-u(2)(i), 
plus n pairs of sm(2)(i) instanton and S'u(2)(2) anti-instanton. For m = 1, different n are 
in different diffeomorphism classes, but in the same homeomorphism class. The standard 
sphere corresponds to n = 0, m = 1, which also corresponds to one instanton of 5^(2) (1) on 
S"^. For n > 0, ??i = 1, it is the exotic sphere, and it also corresponds to n + 1 instantons 
of S'u(2)(i), and n anti-instantons of S'u(2)(2). The case for m = 2, and other general m, are 
also very interesting. 

We propose to realize exotic sphere fibrations (El [3]) by geometric phases in condensed 
matter systems, cold atomic systems or molecular systems. These exotic fibrations can be 
mapped to band structures of those systems with nontrivial band topologies in the systems. 

We start with wavefunction \ip } , where I labels bands in band structure, or labels state 
vector bases, and define gauge field associated with the geometric phase in parameter space, 

/1«^^ = -^(#({U)|5«J^^({U)>, (4) 



A^^^'^^ = -t{ri{U)\d^i^'i{U)) , (5) 

where i,j, a, (3 label different bands, or state vector bases. ^^ are parameters for the wave- 
functions. The i,j labels a double band and takes values 1 or 2. The a, (3 labels a dif- 
ferent double band and takes values 1 or 2. We have written the states in orthonormal 
basis. When the two bands are degenerate or almost degenerate, the geometric phase be- 
comes a non-abelian matrix- valued quantity. ^^ are the coordinates on the parameter space. 
(^15 ^25 ^35 ^4) parametrize four dimensional manifold. Since there are two su{2) gauge fields, 
the band structure is such that there are four occupied bands, and there are two su{2) gauge 
fields associated with the two double bands in the parameter space. The two bands in the 
double band will be degenerate at some points of the parameter space, the ^ space. 

The four-parameter space can be momentum space, or it can be momentum space to- 
gether with extra parameters, or it can be other parameters in the model. It is good to 
embed a three dimensional experimental system into a four dimensional parameter space. 
For example, the four dimensional parameter space (^1,^25^35^4) can be three dimensional 
momentum space {ki, ^2, ks) with an additional parameter ^4. ^4 can be a parameter in the 
model Hamiltonian or the effective model of the experimental system. ^4 can also be k^, or 
u, or a parameter in the model Hamiltonian. It can also be two dimensional momentum 
space {ki,k2) with additional parameters ^4,^3- It can also be real space, together with 
additional parameters. 

The geometric phase is sometimes also called the holonomy. Under the adiabatic evolution 
along a closed path in the parameter space, the state vector will come back to itself up to 
an extra unitary phase factor, given by the integration of path ordered exponential of the 
Berry holonomy along that path. 

We look for experimental quantities that can be realized from the structures of ([2], [3]). 
The field strengthes are 

^i;^'' = -^ {9a^' \9b^') + ' {9b^' \da^') + ^ {^' \9a^') {^' \9b^') - ' {^' \9b^') (^' \da^') > 

(6) 

(7) 
a, 6 label the parameter space, and i,j label the state vector bases or bands. The bracket 
denotes the inner product of wavefunctions in the Hilbert space. One interesting situation 



is when ^^ = ka, a = 1, 2, 3, and ^4 is an additional parameter, and this correspond to the 
momentum space of 3D materials. Summation of the scripts are assumed in the notations. 
Moreover, there are second Chern numbers C2 , C2 of these two gauge fields. 
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^(2) _ 1 
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j (trF(i) A F«) = ^,jd'^ ea.c.(trFi;)Fi])), (8) 

j (trF(^) A F(^)) = ^,jd'i e^UtrFJ^fF^f), (9) 



where the integration is on the parameter space (^1,^25^35^4) &^d c?^^ = d^id^2dizd^A- The 
second Chern numbers are 

C2 = n + m, C2 = —n. (10) 

The non-trivial Chern numbers correspond to the non-trivial topology of the band structure. 
It is an S^ fibration of four-sphere. The su{2) fibration on the 4d parameter space realizes 
a higher dimensional manifold. The ([HI E]) can also be interpreted as four-form magnetic 
monopole fluxes on the parameter space. 

Realization and interpretation of the two su(2)'s and the parameter space can be diverse, 
in condensed matter and atomic systems. This type of fibration are still abstract. It can be 
mapped to band structures in different possible systems. It has relevance to both electron 
band structures and photon band structures. It is possible to have this band structure in 
certain synthesized materials. The Chern numbers may correspond to the number of edge 
states. One of the most interesting situation is m = 1. To exhibit Cg — n -|- 1, C2 — — n, 
the material may have n + 1 right chiral edge states and n left chiral edge states. In the 
context of photonic crystals, it may have n + 1 uni-directional right-moving edge states, and 
n uni-directional left-moving edge states. These topological quantum numbers may manifest 
themselves in terms of the number of edge states, and may contribute to conductivities and 
transport properties of the materials. 

In the situation that the fourth parameter ^4 can be integrated out, the expressions can 
be reduced to Chern-Simons integrals, because of the relation. 
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e.a.cVM^y.A'^i^ + ^,^KA\,A%] = -e^Ut^F^'F'^'). (11) 



Therefore the reduced integration keeps the information of the four dimensional integral. In 
the case of momentum space {ka, kb, kc) plus ^4, the integration after reducing on ^4 is over 
the Brillouin zone. 



The four dimensional topological insulator can be characterized by Z invariant. In 4D it 
may be characterized by C2 G Z. So here the system with Z © Z invariant may be mapped to 
doubled topological insulators in four spatial dimensions. Upon a reduction on the fourth 
parameter ^4, the system becomes a doubled three dimensional topological insulators. The 
three dimensional topological insulator have been characterized by Z2 invariant.-"- The ele- 
ments of Z2 correspond to odd number or even number of Dirac points, which is related to 
the global property of the Brillouin zone.-"- The odd class are topological insulators, and 
the even class are normal insulators. 

The fibration may be related to axion electrodynamics. An effective axion term can be 
induced in several ways^^^i in three spatial dimensions and one time dimension. Similarly, 
in four spatial dimensions, there are effective Chern-Simons action that can be induced^ in 
four spatial dimensions and one time dimension. The Chern-Simons term Ai^e'^^^^^d^AjjdxAp 
is the one-loop effective term arising from integrating the fermion loop and the coefficient is 
given by the second Chern number in the momentum space of the fermion. Up on reduction 
on K = 0:4 direction, A^ becomes 6 field. 

In three spatial dimensions and one time dimension, the systems can effectively have 
axion term, or 6 term, and have magneto-electric effects, 

S = y|^(4'^ + c^?) j dtd\{et^''^^d^A,dxAp), (12) 

^ = (cW+cf))^ = m^, (13) 

where Afj_{x) is real space gauge field. When the system has time-reversal symmetry, the 
time-reversal symmetry and gauge symmetry require that 9 = rriTr, where m is an integer. 
This term is proportional to E ■ B, and thus this will give magneto-electric polarization of 
the material. 

If we consider the interface between two materials with different 6, in which case there 
is 61 for one material extending along 2; < 0, and 62 for the other material extending along 
2; > 0. There is a jump 61 — 62 = A6 across the two sides of the interface. The time- 
reversal symmetry can be broken on this interface.—"^ We can use the integration by parts 
{6)dzAxdtAy = —Ax[{dz6)dtAy] up to total derivatives. Since A^ is coupled to the current 



jx via Axjx, then at the interface we have that the induced current is 

jx = (TxyEy, (14) 

J. = J^ +£\ J^ = -Si^, Ji^^ = -^SEy (15) 

£^^-P^^EA^-in,rnfl^^Ey, (16) 
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The two conductivities have opposite signs. From the expression of the second Chern num- 
bers, the conductance can be expressed as 



d5 = ^ I d'k[{i,' \d,i^') {d.i^' |9,#> + \ {^' \d,^') {r \d.^') {^' I^A^^>]x^\r 



^S = TT^ [d'k[{r |5,^0 {d.iJ^ \dxr) + 1 {r i^.^O i^^ i^^r) ir \dxr)]T^\^e^-^\ 

(18) 



where J d^k = J dkxdkydkz, and where there is summation in the i,j, a, /3 labels of different 
occupied bands. The integration of the second Chern class over / x BZ reduces to the 
integration of two Chern-Simons forms on BZ. 

The fibration is also related to band touching phenomena. There are nontrivial topological 
structures of occupied bands. There are also other unoccupied bands. We can diagonalize 
the Hamiltonian in the system, and it can take the form 

Hik) = E,ik) J2 \h k) (2, k| + E2(k) J2 l«' k) («> k| + 5^ E,ik) I7, k) (7, k| . (19) 

i=l,2 a=l,2 7 

We have that Ei,E2 are energy eigenvalues of the two double bands. Those are occupied 
bands. We have written it in orthonormal basis. E^ correspond to other unoccupied bands. 
The second Chern numbers of the occupied bands are n + m, —n respectively. We propose 
that there may exist band structures with associated geometric phase that realize the exotic 
spheres. 

When we tune the parameters in the model describing the material, the band structures 
can be deformed and changed. As long as there is no band touching that occurs during the 
tuning, the individual Chern number of each band remain unchanged, since it is topologically 
invariant. If when tuning the parameters, the band touching happens, then the Chern 
numbers of two bands that touch can change their individual Chern numbers. In these 



situations, when tuning the parameters, the two bands first touch and then spht. If we view 
the geometric phase gauge field as a fiber bundle over the parameter space, for example the 
three dimensional momentum space plus an additional parameter, then the transferring of 
Chern numbers during band-touching is a topology change of that fiber bundle. The base 
space is the parameter space, and the fiber is the geometric phase gauge field. The Chern 
number of the band may correspond to the number of edge states. 
Consider the effective Hamiltonian of the two double bands, 

^„„ ^ (MSilliiM) ) ,,,, _ (M^O^iMM)) f . .,k;0, (20, 

where we write it in terms of 4 x 4 Hamiltonian iJ(k;,^), and where u -u = 1, and F denotes 
Gamma matrices and C, here denotes a parameter in the model Hamiltonian, for example 
spin-orbit coupling. If the band touching happens at a point near (k; ^), we draw a surface 
S3 enclosing that point. The Chern number transfer between the two bands is therefore 

n = — - / (w \ujduj A duj A doo) . (21) 

This is the Chern number that is transferred between the two bands. 

Starting from the state corresponding to C2 = n + m, C2 = ~''^) we can tune the 
parameters in the model, for example, the spin-orbit coupling, or strain, or magnetic field, 
so as to change the band structure to make the band touching point between the two 
bands happen. The two bands then split by tuning these parameters. When the bands 
touch and then split, they transfer Chern number Ac2 = n, and then the state becomes 
C2 = n' + m, C2 — — n', in which n' = n + n. From the point of view of instantons on the 
base space of the exotic sphere fibration, this transition between the two bands correspond 
to the transfer of h instantons between the two su{2) gauge fields. This also corresponds 
to the topological change transition between exotic sphere fibrations with different n. The 
same procedure of band touching and Chern number transfer can be performed many times. 
This generates different fibrations with C2 = n + m, C2 = —n, and exotic spheres with 
C2 = n + I, C2 = —n. The change of Chern numbers can be mapped to the topology 
change of the fiber bundle. These transitions may relate different n with fixed m. 

Those states with fixed m = 1, fixed n belong to the same diffeomorphism class, and 
they may correspond to the situations that the topology of the band structure does not 
change, and in particular without band touching and Chern number transfer. Those with 
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fixed m = 1, but different n, are in the same fiomeomorpliisni class, but for different n (mod 
28) are not in tlie same diffeomorpfiism class, and tliey may be non-adiabatically connected 
by band toucliings. 

Tliose witli different m are not in the same homeomorphism class and may not be non- 
adiabatically connected. Since the two bands touch and split, the total Chern numbers of 
the two occupied bands, which is m, is topologically invariant and a conserved quantity of 
the total system. There are changes of the individual Chern numbers of each band, and 
there is a Chern number transfer between the two bands. 

Many electron systems can be described by effective Hamiltonian that is quadratic in 
the Fermi field. In some situations there may contain terms that are quartic in the Fermi 
field. The quartic terms can give radiative loop corrections to the quadratic terms. In some 
situations that a mean field theory can be applied, and quartic term may be substituted by 
a quadratic term by a mean field approximation. 

The system may be related to diverse systems, for example, topological insulators, quan- 
tum hall systems, models of Weyl semimetals, semiconductors, photonic crystals, and cold 
atom systems. 

III. THEORETICAL DOUBLED TOPOLOGICAL INSULATOR MODEL 

The 3D topological insulator materials include for example the Bi2Te3, Bi2Se3, 
Sb2Te3, and Bii_a;Sba; (for example Bio.gSbo.i). The alloy Bio.gSbo.i has five Fermi crossing 
points, and its theoretical model is more complicated than the situations of Bi2Te3, Bi2Se3, 
Sb2Te3, which have a single Dirac point. The topological insulator materials can be known 
by measuring whether there is odd or even number of Dirac points, which is related to the 
global property of the Brillouin zone. 

The 3D Topological Insulator materials, for example^"- Bi2Te3, Bi2Se3, Sb2Te3, can be 
described by a simplified effective 4x4 model Hamiltonian, near the level-crossing point. 
They have large bulk gap of order (1 ~ 3)xl0~^ eV. The model with Sel-Bil-Se2-Bil'-Sel' 
crystal structure have been studied in detail.- For example, for a Bi2Se3 crystal, in the 
effective model of the four band model, the basis of the state vector is 

|Pl+,t), \P2-A), \Plt.i). 1^2;, i) (22) 
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where -Pl^, P2~ are two p orbitals in the situations in a Bi2Se3 crystal. Near the level- 
crossing point, the two bands touch each other. In this effective model, the four basis are 
from two orbitals and two spins. Let's abstract the state vector into 

i+,t), i-,t), 1+,;), I-,;). (23) 

We may introduce another su{2) pseudospin space, in which we enlarge the 4x4 model 
Hamiltonian and 4 vector into 8x8 model Hamiltonian and 8 vector 

|+,l,t), |-,l,t), l+,l,;), |-,1,;), l+,2,t), |-,2,t), l+,2,i), |-,2,;). (24) 

The pseudospin refers to the labels 1 and 2 in (!24|) . There are potentially many different 
kinds of 3D topological insulator materials. It is in principle possible in the future to consider 
Interpenetrating Lattices of two topological insulator materials, or Interpenetrating Lattices 
of one topological insulator and one normal insulator material experimentally. Here, we only 
discuss it theoretically. The model Hamiltonian of the 8x8 model can be described as 



H{k) = eo(k,Ol 



8x8" 



-9ik,0 



I4x4 + 



s'\k,on':\k,o 



(25) 

ria (k,^) and ni (k,^) are unit-norm vectors mapped from (k,^) space. It is a 8 x 8 model, 
doubled from 4x4 model.— The model is analogous to doubled topological insulators, and 
have two second Chern numbers 

c? = ^ / ^'fco?ee'^''^'^^n«9fc,n[^)9fc,n«9fc3n«9^r2« =m + n, (26) 

c?^ = ^J d'kdie'^'^'^n^^d,,nf^d,,n^^d,,nfd^n^:^ = -n. (27) 

These integral representations also give the second Chern numbers. 

In that case, the boundary states may compose oim + n right chiral fermion modes, and 
n left chiral fermion modes. The Hamiltonian density of these states in momentum space 
may be expressed as 

H{k)= Y. ^v^4{^■^)i^^+ Y. H-^](-'^-k)^r (28) 

i=l,...,m+n j=l,...,n 

In general this fibration (|8l [9]) may be mapped to m + n right chiral modes, and n left 
chiral modes. The total helicity number is m. 
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It can be realized in Interpenetrating Lattices of two kinds of insulator materials. In 
the context of topological insulators in three spatial dimensions, the situation with the odd 
number of edge states is topologically robust. Since one can perform perturbations to the 
system and a pair of Dirac cones can be coupled and then gapped after re-diagonalization of 
the Hamiltonian. For the odd number of Dirac cones, such perturbations will always leave 
at least one Dirac cone un-gapped. 

For m = 2, it can be realized in Interpenetrating Lattices of two kinds of topological 
insulator materials. For m = 1, it can be realized in Interpenetrating Lattices of one kind 
of topological insulator material, and one kind of normal insulator material. 

Because of the u{l) symmetry, the electric current in the system is exactly conserved 
current. Similar to the discussion in section |lll the electric charge current is 

(1) _ e^ (^i -^2) (1) . . 
J. - ^ 27r "^2 J^v, l^yj 

(2) _ e^ (^i -^2) (2)^ .3.^ 

at the interface between the Interpenetrating Lattices of the two materials with 6^ , 6^ 
respectively, and another material with 62- 

In the context of Interpenetrating Lattices, the enlargement from 4 x 4 to 8 x 8 is due 
to two types of interweaving lattice site L'-^-', L'-^-'. Measurement associated with particular 
sublattice L^^^ or L*^^) selects the corresponding Chern number. 

IV. PHOTONIC CRYSTALS 

We may connect these fibrations to electron band and photon band. The photonic bands 
are parallel and similar to electronic bands. Photonic crystal with bulk band-gap, and 
gapless edge modes are in some aspects similar to topological insulator. It can have several 
bulk band-gaps. It may have special edge states. Since the material has bulk band-gap for 
the photon, it will forbid the bulk transmission of the photons in certain range of frequencies, 
for example UJ2 < 00 < Ui. There can be surface band that are within the band-gap region 
of the bulk bands. 

There are many ways to engineer photon bands in photonic crystals, and there are typ- 
ically many closely-spaced bands. So there are many possibilities to have several Chern 
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numbers. There are Dirac points near band touching points. The photon bands also have 
geometric phases. The 2D photonic crystals can be made by periodic array of cylinders 
of dielectric medium, with lattice structure in x, y directions. The geometric phase of 2D 
photonic crystal can be defined and its first Chern number is given by integration of the field 
strength of the Berry phase gauge field in the 2D momentum space {k^, ky), for example.— 

The 3D photonic crystals (PhC) and 2D photonic crystals (PhC) have a difference that 
the 2D photonic crystals have extra translational symmetry in z direction. The 3D photonic 
crystals can be made by 3D periodic arrays (or lattices) of dielectric spheres, or alternatively 
by 3D periodic arrays (lattices) of air holes in dielectric medium, or by 3D meshes of dielectric 
medium. There can be a limit that the lattice spacing along z direction is much smaller 
than the lattice spacings in x, y directions, and under such limit it cross over to 2D system. 
One can also define a conductivity of edge modes of photons in 3D. 

The band topology can be realized also in three dimensional photonic crystals. The 
geometric phase can also be similarly expressed 



^•,(k)=Im 



u,{k),B-\uj{k))Vk^u,{k) 



V 



u,{k),B-^{u,{k))u,{k)) (u,{k),B-^{co,{k))u,{k) 



(31) 



/ 



AU^) = Im 



(uo,{k),B-\up{k))Vk^U0(k) 



\ 



V 



uM,B-'{ujM)ua{k)) [up{k),B-^{uj^{k))up{k) 



(32) 



/ 



The round bracket denotes contraction of spatial components of the field variables.— Here 
we include non-abelian geometric phases. The i,j labels a double band and takes values 1 or 
2. The a, f3 labels a different double band and takes values 1 or 2. We have assumed almost 
degeneracy coi ~ coj, and coa — i^p- We assume that there are two almost doubly degenerate 
bands in the band structure. 

B~^{r,u) is an 6 X 6 block-diagonal permittivity-permeability tensor 



B~\r,u) = 

/i, 

B^^{uj) has taken into account frequency dependence,— 



eab{r,u) 
f^abir,(^) 



(33) 



B-\r,uj) = B-\r,co) + ujd^B~\r,uj). 



(34) 
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The eab{r,u), /i^f,(r, w) are 3x3 permittivity tensor and permeability tensor, and they 
generally have off-diagonal components. -B^^(a;(k)) is a nontrivial tensor due to the dielectric 
medium, and they have frequency dependence.— 

The ttj(k, r)e^^'^ is the Bloch state of the 6-component complex vector (-Ej(k, r), ifj(k, r)), 
of the electromagnetic fields of the normal mode with momentum vector k and frequency 
Ui(k). 

One can define a Chern-Simons integral, 

/ = ^ / dk,dkydheabc[A^jVk,A'ji + i-AliA^ijA'ji]. (35) 

The integral is in the 3D momentum space {kx,ky,kz) of the 3D photonic crystal, and 
summation of the scripts are assumed in the notations. 
Because of the relation, 

edabcV4At,VtA% + i^AlA\^A%] = ^-eaU^iF'^'F^'') (36) 

the four dimensional integral of 

\ j diJi.d^M, eaU^iF'^'F^') (37) 

can be reduced to three dimensional integral of 

diJi,daA%V,A% + i\aIA\^A''^;\ . (38) 

Therefore the Chern-Simons integral has the information of the four dimensional integral 
that is associated with the C2- 

For particularly engineered 3D photonic crystals, there could be two second Chern num- 
bers C2 ? C2 , whose values may correspond to the numbers of right-moving and left-moving 
boundary states. This type of band structure can also be realized in Interpenetrating Lat- 
tices of two photonic crystal materials PhC and PhC, experimentally. In the latter case, 
C2 ; C2 correspond to the two materials respectively. The Chern numbers may correspond 
to the number of boundary states, or uni-directional one-way propagating states on the 
boundary. This is independent of boson or fermion statistics. Photons can have right cir- 
cular polarization and left circular polarization. In this context, the photon's left or right 
polarization pattern of the boundary photon states would be interesting physical observable. 
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There can be a limit that the lattice spacing along z direction is much smaller than the 
lattice spacings in x, y directions, and the 3D system can crossover to the 2D system. The 
crossover relation between 3D TI and 2D TI in some aspects may be similar to the relation 
between 3D PhC and 2D PhC. 

V. HETEROSTRUCTURE OF TI/NI/TI/NI' 

In this section we discuss heterostructures of periodic units of TI/NI/TI/NI' materials. 
The periodic heterostructure of TI/NI has been devised.— i^^ The TI/NI stands for Topo- 
logical Insulator/Normal Insulator. In periodic TI/NI model,— ^^ there are two kinds of 
interfaces, NI/TI and TI/NI. In the periodic TI/NI/TI/NI' model, there are four kinds of 
surfaces, NI/TI, TI/NI', NI'/TI, TI/NI. These materials are arranged along the z direction 
layer by layer, from top to bottom direction. Here we make the normal insulators NI and 
NI' on the two sides of the same topological insulator to be different. The difference of NI 
and NI' introduces another sm(2) space, the p-space. We make the parameters of two kinds 
of TI/NI, TI/NI' junctions to be different, so in each periodic unit, there are four materials. 
There is experimental method to make the heights of two normal insulators to be the same, 
while making the tunnelings of the surface electrons across the two kinds of normal insula- 
tors to be different. This type of structure can be experimentally performed by many layers 
of periodic heterostructure of thin films. 

In each periodic unit, there are TI, NI, TI, NI' structures. There are surface electrons in 
upper surface and lower surface of the TI materials (in x, y directions). The model Hamilto- 
nians are built from the states on the interfaces, from the surface states of TI's. Adding the 
normal insulator materials can perturb the surface Hamiltonian by adding tunneling terms. 
The periodic structure is along the z direction. The periodic structure then make the states 
become bulk states of the engineered structure. 

The model Hamiltonian is 



-r+(Ai5i,2i + l + A2(5i,2i) + -r (Ai(5i,2i + A25i,2j + ly 



4,i(Cfex2i+l+Cfe^2i). (39) 



This is generalized from the model Hamiltonian of the TI/NI model with two structures 
in each periodic unit.^^ t^^t^.t^ are Pauli matrices, which act on the pseudospin space of 
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upper and lower surfaces, and r"' 



r^ + ir^ , T 



ir^. The k i is the momentum in 



X, y directions. The i and 2j, 2j + 1 label different topological insulator layers. The Hamil- 
tonian (15^ can describe the periodic structure of topological insulators stacked together 
with normal insulators NI and NI' in between, separating the topological insulators. The 
b^a^ + bx<J^ term in (15^ describes spin splitting of the surface states, that can be induced 
by doping each TI layer with magnetic impurities. A^ describes the tunneling between the 
two surfaces of the same topological insulator. Ai and A2 describe the tunneling between 
the surfaces of two nearby topological insulators through the material in the middle, which 
are the NI and NI' respectively. The Ai and A2 parameters have different sizes, and can be 
the same under the limit A2/A1 — )■ 1. The spacing of the periodic structure is d, and the 
total number of the periodic units is A^. The parameters for TI materials are surface Fermi 
velocity vp, bz,bx and tunneling A^. The parameters for NI materials are the tunnelings 
Ai, A2. The bzC^ + b^c^ is a time-reversal symmetry breaking term. 

There are several differences between the configurations here and the configurations in 
previous discussion.—'^ In the configurations there,— "^ there are two structures, the topo- 
logical insulator and normal insulator in each unit. Here, there are four structures in each 
periodic unit, the TI, NI, TI, NI'. Here, we turn on the magnetic term in both z and x 
directions. 

Making a Fourier transformation along the z direction. 









—ikzld 



(40) 



where A^ is the total number of the periodic units, the 3D momentum-space Hamiltonian is 
a 8 X 8 Hamiltonian, 

H = Y,ciH{k)c,, (41) 



Hik) 



e+b.a'+b^a'' 


AJ. 





^^gifc.dj^ 


AJ. 


-e + b^a^'+b^a'' 


All. 








All. 


e+b^a'+b^a'' 


AJ. 


Aae-^'^^'^I. 





AJ. 


— £ + bzO^+b^a 



, (42) 



where e = ^fir(zxcr)-k. The Hamiltonian here is 8 x 8. Because another 5^(2) p-space is 
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introduced when A2 7^ Ai, the Hamiltonian is enlarged from 4 x 4 to 8 x 8. 

We introduce a su{2) space, p-space, where p^, p^, p^ are Pauh matrices. The Hamiltonian 
can be expressed as 

H(k) = [hvpT^ (g) (i X 0-) • k + A,r^ ® !„ + I^ O {b.a' + b^a"")] (g) Ip 

+i[r+Ai + r-A2e''=^>+i ® 1^ + ^[r-A, + r+A^e-^'^^ip-^ ® I., (43) 

where p'^ = p^ + ip^, p~ = p^ — ip^ , and \p is the identity in p-space. 

In the special case, when 6^^^ = 0, making the transformation r^ — t- r^a^, ct^ — t- a^r^, we 
find 

i^(k) = [nt;j.(;z X 0-) ■ k + b,a'] lr®lp + (x'S, (44) 

S = [A,r^' ® I, + ir+(Aip+ + Ase-^'^^^p")^ + ^r-(Aip- + A^e*'^^ V)^]- (45) 

S does not contain a, but only r, p operators, so it commutes with 
[hvp^z X (t) ■'k + bzO'^] lr<S)Ip, and commutes with the Hamiltonian, we can replace it 
by its eigenvalues, T,±. We find 



±S 



± 



± f i(A? + A^) + AlT ^J\{Al - Aiy + (A? + Ai) A2 + 2A2A1A2 cos{Kd) j ,(46) 

and 

H±± = hvF{zxa)-k + b,a' ± E±a'. (47) 

The first ± and second ± subscripts in H±± in P7|) denote the ± in front of S±, and the 
subscript of T,±, respectively. 
The energy eigenvalues are 



E = ±^h^vlkl + {b,±Eiy, (48) 

where ±E-|- is in ( l46l) . 

The Weyl node happens at, for example, when 6^ — S_ = 0. For simplicity, we now denote 
bz as b. The locations are 

kz ='^± ko, (49) 

ko = ^arccos(T^^^=^=^[&^-2(A? + A^ + 2A^)62 + 4(AiA2 + A^)2]i). (50) 
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This state is very similar to a wave-packet of electrons whose kz centers around the particular 
value ^ ± fco- In this situation, the Weyl nodes are centered around {k^, ky, k^) = (0, 0, ^ ± 
ko). So this means that the material only allows the transmission of the state with prescribed 
range of momentum. The limit A2/A1 — )• 1 reduces to previous discussion,—*^ in which case 
both Ai, A2 equal to A^. 

The expansion around the band touching point 6^ — S_ = gives 

H = hvplz X cr) ■ k ± hv^kza"" (51) 

= hvpicr'^ky - a'^k^) ± hv^a^^kz (52) 

where we redefined k^ — {^ — ko) ^ k^, and where 
Near the Weyl points, the dispersion is 



E{k) = ±h^vl{kl + kl) + vlkl (54) 

This is a Weyl fermion in three spatial dimensions, as described in f l52|) . This Weyl fermion 
is half of the components of a Dirac fermion. These Weyl fermions come in pairs, around 
kz =^ ± fco- The two Weyl fermions are separated in the k space and have opposite 
helicities. The heterostructure of topological insulators and normal insulators thus become 
Weyl semimetals. 

If A2/A1 —7- 1, the above reduces to 



S± = ^/ A? + A2 ^ 2A,Ai cos{hzd), (55) 

The branches with different signs can be understood as a phase shift in ^k^d — )■ ^k^d+ir. The 
period when A2 7^ Ai is ^, and when A2/A1 — )• 1 is enhanced to ^ = jj^. When taking 
the A2/A1 -^ 1 limit, Ai = A2 = Aj, 

Al + A^-b^ = 2A,Ai cos{-kod), (56) 

hvo = A,Ai— sin(-A;orf) (57) 

which reduces to previous discussioni^ii^ in the A2/A1 — )■ 1 limit. 
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When we tune the parameters b or the tunnehngs, the band touching occurred and the 
transitions^ between insulators and Weyl semimetals is closely related to the topology change 
of the band structure on the parameter space. 

There is another limit A2/A1 — )■ limit, in which the tunneling across the NI' material 
is taken to zero, so this is the limit when there is TI/NI/TI heterojunction. 

The doping with magnetic impurities can be experimentally performed in, for example,— 
by for example Mn doped Bi2Se3. These heterostructures can be experimentally made by 
using topological insulator thin films, which have been demonstrated in experiments.— 

These models of Weyl semimetals have bulk gapless modes at particular momentum 
vectors, near the above Weyl nodes, and these modes with the particular momentum vectors 
near the Weyl nodes conduct electric current as well as thermal current. These materials 
have electrical conductivity and thermal conductivity which can be measured and useful. 

VI. HETEROSTRUCTURE OF TI/SC/TI/SC 

One can also replace the normal insulator materials discussed in the last section with 
superconductors (SC), and consider TI/SC/TI/SC heterostructures. One can also consider 
the model by changing the NI to superconductor (SC), and construct the model corre- 
sponding to TI/SC/TI/SC heterostructure. The superconductors can introduce couplings 
between opposite spins within the same surfaces (for both upper surface and lower surface), 
|A|e*'^C;[,^c_;^,l+h.c.. The parameters for SC and SC are slightly different. 

Now the Hamiltonian is further enlarged by a K-space. The su{2) K-space is due to 
particle-hole symmetry. The size or dimension of the state vector is doubled due to particle- 
hole symmetry. The periodic TI/SC model has been considered.— The surface Hamiltonian 
of each surface is further perturbed by a superconducting pairing term. 

We again transform the Hamiltonian to momentum space via ( HOl) . We then make the 
transformation r^ — t- r^o"^,cr^ — )■ a^r^. The size of the Hamiltonian is 16 x 16. For 
TI/SC/TI/SC heterostructure, the model Hamiltonian in momentum space is 

^= E 4^^-cfci+ E (i^ie'^4Ai;+h-c.) (58) 

A;;«,J=it,± fc;J=J=±,it 

where the second term is superconducting pairing term with A = lAle**^. Hij is the Hamil- 
tonian without adding the superconducting term, and i, j denote four blocks corresponding 
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to±,± in gZD. 

These four eigenvalues of Hij correspond to sectors of r^, p^ = ±1, ±1, 



H±± = hvpiz X cr) • k + (6^ ± S±)a-' 



(59) 



with the eigenvalues 

{K ± s±) 



b±\^{Al + Al) + AlT J\{A, - A^f + (A? + Al)Al + 2A2A1A2 cos(A;,rf) 



When adding the superconducting term, the 4x4 Hamiltonians are 



(60) 



H, 



i=±± 



^E^m(^' 



vf{z X 0-) • k I, + a%K ± S±)I. + ^(|A|e^^«:+ + \A\e-'^t,-)])^^^^ 



(61) 



J J 



where the basis is t/'^j = {cki-^^Cki]^,c_^^pC_^^^), i = (±, ±). The superscript A in H^^j^ 
denotes including the superconducting term. The first ± and second ± subscripts in H^j^^ 
in (!6T|) denote r^, p^ = ±1, ±1. There are four blocks corresponding to (±, ±). The fe^a^ term 
is a deformation of the quadratic term, and is due to the doping of magnetic impurities in 
the TI. K^, K^, K^ are Pauli matrices, and k+ = /t^' + m^, /t^ = k^ — zk^. The superconducting 
term was from a quartic term, and after a mean field treatment it becomes a quadratic term 
in the Hamiltonian. 

For example, for T^,p^ = —1, —1, 

^-- = iT.'^l-i^M^ X ^) . k 1. + a^[{K - S_)I. + ^(|A|e^^/.+ + \A\e-'^K-)])^^^__. 

k 

(62) 



In the matrix form in K-space, 

[hvp^z X cr) ■ k + (6^ ± S±)cr^ 



Ht=±±(y) 



(63) 



|A|e*^cr^ 
|A|e-*^(T^ [hvF{z X 0-) ■ k + (6^ ± S±)(t^; 

The Hamiltonian in ( l63l) can be diagonalized in k- space. When it is diagonalized in 
K-space, this is given by a transformation 

fe;i,i=(±,±) 
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with the basis of state vector 



wherein 



(^+,fe,^+,{-fc),^-,fe,^-,{-fc))> 



1 1. 



V'+,. = yl^'^^'^C'^^t + 7f ^^^'^^--^4' (65) 



^-,u = -^f ^~'^^^^^t + 7f ^'^^^-M4' (66) 

It can again be diagonahzed in K-space. For t^ , p^ = —1, —1, for k^ = ±1, 

H^^ = hvpiz X 0-) ■ k + a%b, - S_ ± |A|), (67) 

where the superscript ± in H denotes k^ = ±1, and ±|A| in f l67|) correspond to k^ = ±1. 

Similarly for the four sectors t^, p^ = ±1, ±1, 

H^^ = hvpiz X cr) ■ k + a'{b, ± S± ± |A|). (68) 

|A| effectively shifts b^ ± Hj-. The eigenvalue is 



E(k) = ±sjh?vl]^l + {b, ± S± ± |A|)2. (69) 

where 6^ ± S± are in (l60l) . 

The Bogoliubov-Weyl nodes near the band touching points are located at for example 

6^ - S± ± |A| = 0. For 6^ - S_ ± |A| = 0, they are 

A:, = ^ ± kt (70) 

k^ = ^arccos(T^^^=^=^[(&±|A|)^-2(A? + A^ + 2A^)(6±|A|)^ + 4(AiA2 + A^)^]^) 

(71) 

and the Fermi velocity is 

_ AgA,A4sin(fcod) 

^"°-A^-(6±|A|)2 + i(Af + Ai)- ^^'^ 

We have thus found the energy eigenvalues, the locations of Weyl nodes, and the Fermi 
velocity near the node. 

Interestingly, Weyl superconducting phases can also be realized by triplet pairing 
phases.— li^ 
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VII. OTHER RELEVANT MATERIALS AND DISCUSSION 

The su(2) geometric phase with two Chern numbers may occur also in other condensed 
matter systems or atomic and molecular systems. There are related discussions in other 
possible situations.—"— The Semiconductor model includes at least light hole (LH) and heavy 
hole (HH) bands. If both LH and HH are doubly degenerate, they can also exhibit su{2) 
geometric phases with two Chern numbers, in the d-space. There are related discussions^^'— 
to these aspects. It is possible to exhibit general value of Chern numbers, by for example 
multilayer heterostructures of particular semiconductors. 

The cold atomic system with cold atom of a larger total angular momentum oi F = 3/2, 
can also exhibit su{2) geometric phase in the parameter space of paring condensates, for 
example.— It may also be relevant to Interpenetrating Lattices of two optical lattices of cold 
atom systems. 

It may be interesting to see whether these fibrations can be realized concretely in these 
experimental settings. 
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